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ABSTRACT 

We present an revised table of 390 Galactic radio supernova remnants (SNRs) and their basic pa- 
rameters, Statistical analyses are performed on SNR diameters, ages, spectral indices, Galactic heights 
and spherical symmetries. Furthermore, the accuracy of distances estimated using the X-D relation 
is examined. The arithmetic mean of the Galactic SNR diameters is 30.5 pe with standard error 1.7 
pe and standard deviation 25.4 pe. The geometric mean and geometric standard deviation factor 
of Galactic SNR diameters is 21.9 pe and 2.4, respectively. We estimate ages of 97 SNRs and find 
that the supernova (SN) birth rate to be lower than, but within 20 of currently accepted values for 
SN birth rate. The mean spectral index of shell-type SNRs is —0.51 + 0.01 and no correlations are 
found between spectral indices and the SNR parameters of molecular cloud (MC) association, SN 
type, diameter, Galactic height and surface brightness. The Galactic height distribution of SNRs is 
best described by an exponential distribution with a scale height of 48+4 pe. The spherical symmetry 
measured by the ovality of radio SNRs is not correlated to any other SNR parameters considered here 


or to explosion type. 
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1. INTRODUCTION 


For their role in Galactic and interstellar medium 
(ISM) evolution, the study of supernovae (SNe) and 
their remnants (SNRs) is vital. In order to better un- 
derstand SNRs, their basic parameters such as distance, 
age, spectral index ete., as well as an accurate count are 
necessary. 

Green (2019) gives the current total of Galactic SNRs 
as 204. While the predicted total number of Galactic 
SNRs is a few thousand (Supan et al. (2018), Ranas- 
inghe & Leahy (2022), observations and identification 
of SNRs are greatly affected by selection effects. 

Several hurdles lead to a smaller sample of Galactic 
SNRs (< 294) with which to perform statistical stud- 
ies. Out of the observed SNRs, many do not have а 
distance or an age because of the difficulty in estimat- 
ing them. The distances to SNRs are estimated using HI 
absorption spectra and maser and molecular cloud (MC) 
associations. When distances using these methods are 
unavailable, distances are often estimated using the -D 
relation (surface brightness ()- physical diameter (D) 
relation). The main drawback of the Y-D relation is the 
large scatter in the X-D plane leading to distances that 
vary from the true distance by about an order of magni- 
tude (Green 1991, 2005). Xu et al. (2005) presented five 


methods to estimate the ages of SNRs. Only a subset of 
SNRs have parameters needed to estimate ages leading 
to even smaller samples. 

While the Galactic SNR sample is far from complete, 
many authors have performed a statistical study on 
them (e.g., Green (1991); Xu et al. (2005)) where the 
the sample sizes range from 174 to 234. The sample 
we utilize for this work is an up-to-date list of Galac- 
tic SNRs which includes uncertain or newly identified 
SNRs. The 390 sources we use in this work is at least a. 
~ 60% increase from previous studies 

In Section 2 we present the updated catalogue of 390 
SNRs SNR with updated distances and other parame- 
ters, In Section 3, the analysis and discussion on the 
radio size distribution of Galactic SNRs, SNR distances 
and У-Р relation, SNR ages, SNR spectral index dis- 
tribution, galactic height distribution of SNRs and SNR 
spherical symmetry is given. Finally, the conclusions are 
given in Section 4 


2. SUPERNOVA REMNANT SOURCE LIST AND 
THEIR PARAMETERS 
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For the source list, we use the catalogues presented by 
Green (2019)! and Ferrand & Safi-Harh (2012)?. Addi 
tionally, an independent literature search was performed 
to search for newly identified SNRs. The catalogue pre- 
sented by Green (2019) has 294 SNRs *. The catalogue 
presented by Ferrand & Safi-Harb (2012) has 383 SNRs 
which includes the 204 SNRs from the Green (2019) 
catalogue. The SNR list presented by Ferrand & Sali- 
Harb (2012) includes new and uncertain SNRs. The 
independent literature search produced 7 newly identi- 
fied SNRs bringing the total to 390 sources, However, 
for this study we do not include the SNR candidates 
(c.g. presented by Anderson et al. (2017)) as a caution 
since it may add more contamination than information 
about real SNRs. The Galactic SNR list and their pa- 
rameters are given in Table 1 (the table in its entirety is 
in the online version of the manuscript). The references 
for the spectral indices, MC interactions, ages and SN 
types are in the online version of Table 1. The source 
numbers (Table 1) first follow the Green (2019) cata- 
logue (1 — 294) followed by the remaining Ferrand & 
Safi-Harb (2012) catalogue (295 — 383) SNRs and finally 
the newly identified SNRs (384 — 390). 

The angular sizes, the SNR types and the 1-GHz flux 
densities are from the catalogues. The radio spectral 
indices were found in the catalogues as well, however, 
performing a literature review we verified the spectral 
indices. For some cases where the 1-GHz flux densities 
were not given, we infer them from the spectral indices 
and flux densities at other frequencies, 

Of the 390 SNRs, 270 SNRs have a radio spectral 
index while 110 do not. Here the radio spectral index 
(o) is defined as S, ox v^, where S, is the flux density 
at frequency v. However, the radio spectrum for some 
SNRs is not well described by a simple power-law, ei- 
ther because the radio spectrum is described by multiple 
power-laws or because of the variations of the spectral 
index across the face of the SNR (Green 2019). There 
are 7 such SNRs in the sample where we denote the 
varying spectral index as “varies”, Furthermore, there 
are 3 SNRs with spectral indices but the values are 
uncertain (denoted with a ‘?'), The 1-GHz flux densi 
ties are given for 183 SNRs and uncertain flux densities 
(denoted with а ‘?’) in Table 1 are given for 108 SNRs 
(none for 99 SNRs), bringing the total number of SNRs 
with flux densities to 291 

The distances for this study are literature estima- 


1 isttp://www.mrao.cam.ac.uk/surveys/smrs/ 
2 nttp://www.physies.umanitoba.ca/snr/SNReat 


Some of the SNRs in the catalogue have been removed (e.g. 
916.8 — 1.1) as more information became available. 


tions (see Ranasinghe & Leahy (2022); the methods of 
the distance estimations and references therein). We 
have added 6 more SNRs (uncertain) to the list with 
distances bringing the total number of sources with 
distances to 221. Green (2019) removed 5 SNRs from 
the catalogue because Anderson et al. (2017) presented 
evidence that they may have been confused with HII re- 
gions. We have included these uncertain SNRs, namely 
620.4 + 0.1, G21.5 — 0.1, G23.6 + 0.3, G59.8 + 1.2 and 
G65.8 — 0.5 in our source list. It should be noted that 
four of these SNRs have distances. The other 3 sources 
with distances that have been included in the source list 
are G190.2+ 1.1, 6192.8 — 1.1 and G359.9+0.0. For the 
SNR (uncertain) G359.9 + 0.0, its association with Sgr 
А* places the SNR at the Galactic centre (GC) (Wang 
et al. 2006). 

‘The ages of SNRs presented in Table 1 are literara- 
ture values except for the ones that are followed by а + 
The ages of SNRs with a ‘*' were estimated using the 
software presented by Leahy & Williams (2017) which 
uses X-ray emission to obtain more accurate age esti- 
mates. There are 115 SNRs with well constrained ages 
and 6 SNR where the ages are given as a lower or upper 
limit. For this work we estimated ages for 95 SNRs. 
А description of the method and parameters used to 
estimate the ages are given in Section 3.3. 


3. ANALYSIS AND DISCUSSION 


The statistical analysis and discussion of the results 
for the 390 SNRs listed in Table 1 are presented here, 
First, we explore the size distribution of the Galactic 
SNRs and the accuracy of the X-D Relation for the 
Galaxy. Next the age and spectral index distributions 
are investigated. Finally, the Galactic height distribu- 
tion and spherical symmetry of SNRs are investigated 
to ascertain whether trends exist for a particular class 
of SNRs 


3.1. The Radio Size Distribution of Galactic SNRs 


There are 221 Galactic SNRs (including uncertain 
SNRs) with distances. Out of the 221, 6 SNRs gives 
the distances as a lower or an upper limit. We exclude 
GLA — 0.1 since its location near the GC lacks strong 
evidence. The radii of the remaining 214 SNRs are given 
in Table 1. In case where the angular size of a SNR is 
given as the major and minor axes, the average radius 
(rag) was estimated using rag = утау ттт, Where 
Tmay and tmin are the semi-major and semi-minor axes 
measured in radio, respectively. For the error of the 
mean radius, we adopt the error estimated using semi 
major axis, The probability density function (PDF) for 
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Figure 1. Left: The smoothed probability density function of the SNR diameters. The red solid line is arithmetic mean 
(30.5 + 1.7 pe) and the vertical black dashed line is the geomatric mean of the diameter (21.9 pc) of the Galactic SNRs. The 
histogram is in the background (in green). Right: The cumulative distribution function of SNR diameters of galaxies. The 
black solid line- SNR diameters from this work (sample size- 214), The red dashed line- M31 from Lee & Lee (2014) (sample 
ize- 156), The blue dash-dot line- M33 from Long et al. (2010) (sample size- 137), green dotted line- M83 from (Dopita et al. 
10) (sample size- 47), magenta solid line- LMC from Bozzetto et al. (2017) (sample size- 59) and cyan dotted line- SMC from 
Filipovié et al. (2005) (sample size- 21). 


the radio diameters is shown in Figure 1 (left panel). diameters are skewed with a bias towards SNRs with 

‘The statistics on the Galactic SNR diameters are given smaller physical diameters. The geometric mean of 
in table 2 and as a comparison, statistics on SNR di- Galactic SNR diameters is comparable to the M83 SNR 
ameters from other galaxies are included as well. The diameter means, 


arithmetic mean of the Galactic SNR radio diameters is 
30.7 + 1.7 pc. The estimated error here is the standard 


error, given by a. / V/ri, where ø is the sample standard ‘The SNR samples of most of the galaxies in the above 
deviation and n is the number of elements in the sam- studies have a striking similarity, which is the lack of 
ple. The median of the sample is 24.3 pe. The geometric SNRs with small physical diameters. Only a few SNRs 
mean of Galactic SNR radio diameters is 21.9 pe. The (<~ 5) in each sample have SNRs where the diameter is 
geometric standard deviation factor (og = exp™'"(P)) < 10 pe. Past authors (Long et al. 2010; Badenes et al. 
of the SNR diameter sample is 2.4. 2010) have stated that this deficit is likely a real one. 

For M33, the mean and median SNR diameter are There are ^ 35 SNRs with diameters < 10 pe in our 
~ 50 and 44 pe, respectively (Long et al. 2010) (see sample. While the relatively large number of SNRs with 
Table 2). Lee & Lee (2014) gives the mean diameter small physical diameters are present in the Galaxy, the 
of SNRs/candidates as ~ 48 pe found in M31. The sample of SNRs/candidates presented by Winkler et al. 
mean diameter of confirmed SNRs in the Large Magel- (2021) (and references therein) for M83 has almost a 100 
Janie Cloud (LMC) is 39 + 4 pe (Bozzetto et al. 2017) SNRs with small physical diameters that are likely ex- 
and the Small Magellanic Cloud (SMC) SNR/candidate panding in locally dense regions of the ISM. Most of the 
mean diameter is 40.8 pc (Filipovié et al. 2005). While SNRs in other galaxies have been identified using Ha 
the mean SNR diameters from M31 and M33 are con- and [S Ш lines, where the identification of SNRs with 


sistent with each other, they are larger than that of the smaller diameters is not expected (Winkler et al. 2021) 
Milkyway SNR mean diameter. Similarly, the mean di- Galactic SNRs are generally identified using different 
ameters of LMC and SMS are consistent with each other criteria (see Brogan et al. (2006)) and are observed 
and larger than the Galactic mean SNR diameters. The at higher spatial resolution, which explains the higher 
mean diameter of SNRs in M83 is 22.7 pe (Dopita et al number of Galactic SNRs with smaller diameters. 
2010) and is smaller than the Milkyway mean SNR di- ‘The cumulative distribution function (CDF) of Galac- 
ameter, tic SNR diameters is shown in Figure 1 (right panel) 
The geometric mean of the Galactic SNR diameters For comparison, the CDFs of five other galaxies (M31, 
(21.9 pc) is smaller than that of the geometric mean M33, M83, LMC and SMC) are given as well. The 
of SNR diameters from galaxies M31, M33, LMC and Kolmogorov-Smirnov (KS) test was performed on the 


SMC. As seen in Figure 1 left panel, the Galactic SNR diameter distributions and showed that the Galactic 
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Table 2. Statisties on SNR diameters of galaxies 


Galaxy — Sample References Arithmatie Standard Standard Median Geometrie Standard 
Size Mean (pc) Dev (pe) Error” (pe) (pe) Mean (ре) Dev Factor 

Milkyway — 214 Table 1 30. 254 LT 21.9 

M31 156 Lee К Lee (2014) 48. 19.0 15 A 

M33 137 — Longetal (2010) — 49.6 28.8 25 42, 

M83 АТ Dopita et al. (2010) 22.7 10.3 15 20.5 

LMC 59 Bozzetto et al. (2017) 39.0 21.3 28 31.6 

SMC. 21 Filipović et al. (2005) 40.8 18.4 40 375 


Nore—a - Definitions are given in Section 3.1 


SNR diameter distribution is not consistent with any of 
the other distributions. The discrepancies between the 
diameter distributions could be real but are affected by 
the different methods of SNR identification. 


3.2. Supernova Remnant Distances and X-D Relation 

While kinematic distances to many SNRs have been 
estimated using HI data, MC or maser associations, op- 
tical extinction, ete., there are a considerable number 
of SNRs with no distances (~ 1/3 of the Green (2019) 
catalogue and ~ 1/2 of the Ferrand & Safi-Harb (2012) 
catalogue). When there are no other means to estimate 
the distances to SNRs, the surface brightness (X)- phys- 
ical diameter (D) relation (E-D) has been employed. 
The main idea is that the surface brightness of a SNR 
changes as the SNR expands or as the diameter increases 
with time. The X-D relationship between the radio sur- 
face brightness at a frequency v and the diameter of a 
SNR is given by 


E,- AD^. a) 


With known SNR distances (from other methods) as 
calibrators, an empirical X — D relation is determined 
to estimate the unknown SNR distances. It should be 
noted that this relation is mainly applicable for shell- 
type SNRs. However, the majority of the Galactic SNRs 
are shell-type SNRs (~ 80% of the Green’s SNRs). 

Many authors have used the X-D relation in the past 
io estimate distances to SNRs. Both Case & Bhat- 
tacharya (1998) & Pavlovié et al. (2014) have presented 
power-law index values (8 = 2.4 and ~ 5.2, respectively) 
and new distances. Their calibrator sample sizes were 
36 and 65 SNRs, respectivel 

Out of the 214 SNRs with diameter estimations, 29 do 
not have 1-GHz flux densities. Out of the remaining 187 


SNRs 139 SNRs are shell-type. For our sample, we have 
included the uncertain shell-type SNRs as well. Figure 3 
shows the distribution of the SNRs in the X-D plane and 
а clear decrease of У with increasing D, on aver- 
age. However, there is a large scatter in the distribution 
and for some SNR diameters there is a wide range of sur- 
face brightnesses for different SNRs. For an example, at 
а ~ 80 pe diameter, the SNR surface brightness has a 
range of values spanning 3 orders of magnitudes. This 
scatter in the -D plane has been pointed out by Green 
(1991) and as recently by Pavlović et al. (2014). Not all 
SNRs evolve the same and is affected by parameters such 
as explosion energy, ambient density, ete. Therefore, a 
large scatters in the X-D plane are expected (Pavlovié 
et al. 2018). Furthermore, the importance of ISM den- 
for the strength of radio luminosity of SNRs was 
found by Leahy et al. (2022), which is not predicted in 
analytic models for radio emission, While distance es- 
timations using a E-D might be more accurate in other 
galaxies (e.g. Bozzetto et al. 2017), it is unlikely a single 
relation could be valid for our Galaxy. 

As Green (1991) stated, it may be possible to improve 
the X-D relation if the sample was taken from a par- 
ticular class of SNRs (e.g. SNRs associated with MC) 
We have investigated the relation taking samples by MC 
associations, SN type and Galactic height (2) and found 
none of the X-D distributions provided a better relation. 
‘The power law index, 3 from Equation 1 as calculated 
in SNR radio emission models, explicitly depends on the 
spectral index, а (Pavlovié et al. 2014). We investigated 
the У-Р relation for different subsets of spectral indices 
and found that the the different sets of fixed spectral 
index ranges did not show any reduced scatter with re- 
spect to that for the whole set. This may be due to the 
fact that the spectral indices are an average across the 
face of SNRs. Regardless, the samples for each spec- 
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Figure 2. A comparison between the literature distances and distances obtained using the X-D relation. The red circles denote 


the literature distances to SNRs from Ranasinghe & Leahy ( 
Case & Bhattachas 
in the x-axis labels corresponds to SNR numbers from Table 1 


tral index (40.02) are too small (~ 10 or less) to make 
а definite conclusion, The X-D distribution in Figure 3 
shows a clear downward trend but only for young SNRs. 
‘The older SNRs that are > 20 kyr in age (diameters be- 
tween 100 & 200 pe) are scattered showing no discernible 
trend. The ages of SNRs come from literature and esti. 
mations from this work (see Section 3,3 for details) 
The У-Р relation has been used as a first guess when 
there are no other methods of distance estimations are 


s could differ from. 


available. However, the X-D distan 
the actual distance by an order of a magnitude and the 
large scatter in the Galactic X-D distribution leads to 
large errors. Figure 2 shows the literature distances from 
Ranasinghe & Leahy (2022) and B-D distances (Case & 
Bhattacharya 1998 & Pavlović et al. 2014). As seen 
in Figure 2, only a few X-D distances agree with the 
literature values. For this comparison we only use the 


2). The distances obtained from the X-D relation presented by 
(1998) and Pavlović et al. (2014) are denoted by blue ‘x's and green diamonds, respectively. The numbers 


Green (2019) SNRs (sample size: 204). Figure 4 left 
panel shows a Venn diagram with the number of SNRs 
in each study, The right panel of Figure 4 shows the 
percent difference between the X-D and literature dis- 
tances. There are only 20 distances that are within 5% 
of each other. The majority differ at least by 20% (on 
average differ by 50%). Thus, while X-D distances may 
give an idea of SNRs’ location, the reliability of the re- 
sult is questionable, 


3.3. Supernova Remnant Ages 

The five main methods of estimating SNR ages use i) 
age of associated pulsar ii) the explosion energy of the 
progenitor and the energy loss rate of the SNR, iii) the 
expansion velocity and radius of the SNR, iv) the ther- 
mal temperature measured from X-ray observations, and 
v) the remnant’s measured spectral break due to syn- 
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chrotron losses and its magnetic field (Xu et al. 2005). 
However, the uncertainties of the ages are large (up to 
an order of magnitude) and the parameters needed to 
estimate the ages for many SNRs are unknown. 

Out. of the 390 SNRs in our sample, there are 141 
SNRs with literature ages, where 9 are given as a lower 
or upper limit. Of the 141 SNR, Sedov ages to 22 SNRs 
were presented with assumed parameters (Ey and по). 
The remaining 119 literature ages were estimated using 
more accurate methods (e.g. X-ray data, shock veloc- 
ity). 

We estimate ages for 75 SNRs with no previously pub- 
lished ages and present revised ages to the 22 SNRs with 
previously published Sedov ages (97 SNR ages in total). 
To estimate the age, we assume that the shock radius 
of each SNR to be the average radio radius. With no 
additional information (e.g. distance, radius), 174 SNRs 
in the sample have no ages. 

To estimate the ages we use the software presented by 
Leahy & Williams (2017). Here we assume the explo- 
sion energy is assumed as, Ey = 0.75 x 10° erg, a value 
between the standard explosion energy of Ey = 10°! erg 
and the logarithmic mean explosion energy of 0.5 x 1071 
erg presented by Leahy (2017) (for LMC SNRs). If the 
SN explosion type of the resulting remnant is known, 
the ejected mass was assumed to be 1.4Mo for type Ia 
SN event and 5Mg for core-collapse (CC). The majority 
(about 85%; Tammann et al. 1994) of the SNe are CC, 
where the progenitor is a young massive star. Therefore, 
if the progenitor type of a SNR is unknown, we assume 
the ejected mass to be 5Ма. Furthermore, we assume 


the interstellar (ISM) number density, по to be 10 em~* 
if the SNR is associated with MC and 0.01 cm if it 
is not. For probable or unknown MC associations, we 
take the по to be 1 cm^?. However, when available, we 
use the literature number density values. The ages were 
estimated for an ejecta power-law index, n = 7 and the 
circumstellar medium (CSM) power-law index, s = 0. 

‘The assumed explosion energies and ISM densities 
could differ by a few orders of magnitudes from the ac- 
tual values. Leahy & Ranasinghe (2018) and Leahy et al 
(2020) presented ages for 58 SNRs, estimated using X- 
ray shock temperatures and emission measures and pro- 
vided explosion energies and ISM densities. As a test to 
evaluate the accuracy of SNR age estimations described 
above, we re-calculated the ages for these 58 SNRs using 
the fixed assumed explosion energies and densities (Eq 
= 0.75 x 10°! erg and по to be 1 cm^? ). Most of the 
ages differ by factors of ~ 1/4 to ~ 4 times (on average 
1) compared to the ages estimated using X-ray data. In 
some cases, the SNR ages differed by an order of mag- 
nitude. However, considering that majority of the SNR 
ages fall between +50%, we adopt it as the uncertainty 
of SNR ages estimated in this work. 

It should be noted that some of the ages estimated 
using the above method are likely incorrect because of 
incorrect assumptions. For an example, G65.8 — 0.5 at 
а distance of 2.4403 kpe gives an age of 500 yr for an 
Eo = 0.75 х 10" erg and по = 1 cm^?. The age for 
G65.8 — 0.5 is likely > 500 yr with either a lower ex- 
plosion explosion energy or higher ISM density. If the 
explosion energy is lower by 2 orders of magnitude or the 
ISM density is higher by 2 orders of magnitude, the age 
of G65.8 — 0.5 is 6000 and 3500 yr, respectively. How- 
ever, without additional information, it is difficult to 
constrain the ages for SNRs such as G65.8 —0.5. There- 
fore, we do not discard these ages from further anal 

Figure 5 shows the plot of SNR diameter against their 
age. To determine the diameter-age (D — t) relation we 
have only included the literature ages estimated using 
more reliable methods (e.g Pulsar associations, expand- 
ing velocity, X-ray observation, etc.). The sample con- 
sists of 83 shell-type SNRs that have literature ages. For 
comparison, in Figure 5 we have inchided lines for the 
main stages of SNR evolution, The ejecta-dominated 
stage (green dashed line) has г œ t4/ (Chevalier 1982), 
adiabatic stage (red dotted line) has r œ 2/5 (Sedov 
1959) and the radiative stage (magenta dash-dotted line) 
has r x 22/7 (Cioffi et al. 1988). Here r is the shock ra- 
dius and t is the age of the SNR. The best-fit line for 
the data is given by 
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Figure 4. Left panel: Venn diagram showing the number of SNRs with distan 
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Figure 5. Galactic SNR diameters and their age distribu- 
tion for a sample of 80 shell-type SNRs. The best fit line 
(black curve) is D = (9.52 + 0.92) 19-59?! pe. The green 
dashed line is r x t*/? ejecta-dominated stage, the red dot- 
ted line is г x t?/? adiabatic stage magenta dash-dotted line 
r x A radiative stage. 


Similar to the YX-D distribution, the X-t (surface 
brightness - age) distribution shows a clear downward 
trend, However, a single relation does not fit the data 
because of the scatter spanning 2-3 orders of magni- 
tudes, Nevertheless, the downward trend does confirm 
that the SNR surface brightness diminishes with age on 
average. 

Figure 6 left panel shows the histogram of the lit- 
erature SNR ages in red and all SNR ages (includes 
literature ages and estimations from this work) in our 
sample in grey. Most of the literature ages (> 95%) 
are < 65 kyr. The main reason for the lack of older 
SNRs is that they are not bright enough to have an 


easily measured X-ray spectrum (Leahy & Ranasinghe 
2018). The sparsity of older SNRs (>~ 50) is appar- 
ent in both samples (literature ages and all) reiterating 
the incompleteness of the sample of SNRs with ages. 

Figure 6 right panel shows the cumulative age dis- 
tribution of Galactic SNRs 
order). There are 204 SNRs with ages in the sample 
(excluding ages given as upper and lower limits). The 
best-fit line for the youngest 180 SNR ages gives a SN 
rate of 1 SN per 164 yr with a correlation coefficient of 
~ 0.87. Owing to the incompleteness of the sample of 
SNRs with known ages, the best-fit line of 1 SN per 164 
уг does not describe the trend well (Figure 6 red dotted 
line compared to the data in green). 

‘The youngest 125 (< 14 kyr) and 31 (< 2 kyr) SNR 
ages in the sample give better best-fit lines of 1 SN per 
87 yr and 1 SN per 58 yr, respectively, both with corre- 
lation coefficients of > 0.98. While the best-fit line of 1 
SN per 58 yr is smaller than the lower limit of the SN 
rate presented by Tammann et al. (1994), it is within 
2o. Furthermore, 1 SN per 58 yr rate is consistent with 
the CC SN rate of 1 per 6171 yr presented by Rozwad- 
owska et al. (2021). From Figure б it is seen that the 
incompleteness of the sample for older SNRs increases, 
‘The sample is likely nearly complete for the SNR ages 
< 2 kyr, but the incompleteness is a factor of ~ 2 for 
SNR ages < 14 kyr. 


3.4. Radio Spectral Indices of Supernova Remnants 

A non-thermal spectral index (о < 0) is one the cri- 
teria in identifying a SNR (Brogan et al, 2006). Figure 
7 left panel shows the PDF of the Galactic SNR spectral 
index distribution, We find the mean spectral index of 
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Figure 6. Left: Histogram of literature SNR ages (red) and all SNR ages in our sample (grey). Right: The cumulative 
distribution of SNR Ages (180 youngest sources). The best fit lines: 1 SN in 164 yr (for ages < 50 kyr, red dotted line, sample 
size: 180), 1 SN in 87 yr (for ages < 14 kyr, black dash-dot line, sample size: 125) and 1 SN in 58 yr (for ages < 2 kyr, blue 
dashed line, sample size: 31). The orange line is the SN rate of 1 SN in 40 + 10 yr presented by Tammann et al. (1994) with 
uncertainty shown by the orange shaded region. 
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Figure 7. Left panel: The PDF of the spectral index (a) distribution. SNRs (all types: 270 SNRs in sample) in red (histogram 
and PDF) and shell-type SNRs (208 sample size) in blue (histogram and PDF). The vertical lines show the mean spectral index 
Tor each SNR type. Right panel: The spectral index vs age (175 SNRs in sample). The horizontal lines show the mean spectral 
index for each SNR type (same as left panel). 


all Galactic SNRs to be а = —0.49 + 0.01 with a 0.16 The mean composite SNR spectral index is a = 
Standard deviation. The estimated error is the standard —0.41 + 0.03, which is slightly flatter than the mean 
error. The statistics on spectral indices of SNRs are shell-type SNR spectral index. There are 37 composite 
given in Table 3. The mean spectral index of the shell- SNRs in the sample including 17 uncertain composite 
type SNRs Ва = —0.51 + 0.01. Both mean spectral type SNRs (denoted with a *?)). The composite spectral 
indices are consistent with the theoretical SNR spectral indices range from —0.68 to —0.1 with a standard devia- 
index of —0.5. The mean spectral index of the LMC tion of 0.15. However, it should be noted that the indices 
SNRs is а = —0.52 (Bozzetto et al. 2017) and of the < —0.2 are classified as uncertain composite SNRs. 
SMC SNRs Ва = —0.63 (Filipović et al. 2005). We 

find the mean Galactic SNR (all and shell-typo) spec- The mean filled-center SNR spectral index is a = 
tral indices are consistent with the LMC mean index, —0.34 + 0.08. The spectral index of filled-center SNRs 
however, the SMC SNR/candidate mean spectral index are characterized as being flat where a ~ 0 — —03 


is steeper than the Galactic mean index. (Weiler & Panagia 1978). There are 11 filled-center 
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‘Table 3. Statistics on SNR radio spectral indices 


SNR Mean Standard Standard Median 
type dev error 
АШ -049 — 016 001 — -050 
Shell -051 0м 001 -0.50 
Filled-center -0.34 — 027 008 0.26 
Composite -0.41 015 003 — -039 
Unknown -0.49 — 017 005 — -045 


SNRs and majority of the spectral indices fall in the 
range of a ~ 0 — —0.3. In fact only 3 filled-center SNRs 
have steeper spectral indices (G65.7 + 1.2, G76.9 + 1.0 
and G141.2 + 5.0). If the three steeper spectral in- 
dices are excluded from the filled-center SNR sample, 
the mean spectral index is —0.20--0.03 with a standard 
deviation of 0.10. While the filled-center SNR sample 
size is too small to form a definite conclusion, both mean 
spectral indices (total sample and sample excluding the 
three steeper а) are consistent with the generally ob- 
served filled-center SNR spectral index. 

Additionally, we have analysed the spectral index dis- 
tribution of known pulsar wind nebulae (PWNe). The 
PWNe with radio counterparts were taken from the cat- 
alogue presented by Roberts (2004)*. Similar to filled- 
center SNRs, the spectral indices of PWNe are expected 
to be relatively flat. There are 34 PWNe in the Roberts 
(2004) catalogue that have associated SNRs, where 31 
have spectral indices (3 have varying as across the 
SNRs). The PWNe spectral indices range from —0.68 to 
—0.06. The mean spectral index of the PWN sample is 
—0.36+0.03 with a standard deviation of 0.17. However, 
unlike the filled-center SNR spectral index distribution, 
the PWNe spectral indices appear to follow a normal 
distribution. 

We have investigated the spectral index distribution 
against MC associations, SN type, diameter and sur- 
face brightness and found no correlation between them 
confirming the conclusion formed by Clark & Caswell 
(1976). The age-spectral index plot (Figure 7, right 
panel) shows significant scatter and no apparent trend. 
Younger SNRs are expected to have a flatter spectral 
index (Dubner & Giacani 2015), but the majority of the 
young SNRs (< 1 kyr) in our sample are seen to have 
a steeper spectral index (average of ~ 0.7). The older 


4 https: //www.physies.mnegillea/~pulsar/pwneat. html 


SNRs (> 100 kyr) show an average spectral index close 
to the theoretical value of ~ 0.5. This correlation, while 
weak, seems to be consistent with the LMC SNR spec- 
tral index distribution (see Bozzetto et al. 2017). 


3.5. Galactic Height Distribution of SNRs 


‘The majority of the SNRs are located in the Galactic 
plane and most of them are CC that can be associated 
with star forming regions. We have investigated whether 
there is a correlation between the SNR diameter and the 
Galactic height, z. There is a large scatter in the distri- 
bution and the correlation of size with Galactic height is 
poor (Figure 8). However, while the correlation is poor, 
there is an upward trend, indicating that as the galactic 
height increases the SNR diameter increases. Further- 
more, the scatter in the D — plane is not further re- 
duced when considering only SNRs with or without MC 
associations 


100. 
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Figure 8. The Galactic radio SNR diameters, D vs. Galac- 
tic height, z for the sample of 214 radio SNRs. 


Figure 9 shows the Galactic height distribution of 
SNRs. We fit exponential and normal distributions 
to the data ). The best-fit line 
for the exponential distribution is given by № = 38.8 
exp(—|2|/48.1) and for the normal distribution is given 
by N = 23.9 ехр(-0.5(2/57.3)?). The minimized x? 
values for the 2 models are 28.2 and 41.5, respectively. 
With a degree of freedom of 59, the p-value for the 
exponential distribution is > 0.99 and for the normal 
distribution is 0.96. Therefore, examining the curves 
and p values, we conclude that the exponential distri- 
bution describes the data best. 

For the Galactic height exponential distribution of 
SNRs, most of the SNRs are located near the galactic 
plane. In fact ~ 90% of the SNRs are located at a 
Galactic height < 200 pe. This is consistent with the 
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Figure 9. The Galactic height distribution (z) distribution 
of SNRs (sample size of 214). The best-fit exponential dis- 
tribution, № = 38.8 exp(—|2|/48.1) is the black solid curve. 
‘The blue curve is the best-fit normal distribution, N = 23.9 
exp(-0.5(z/57.3)*). 


massive star and star forming region distribution of the 
Galaxy (Urquhart et al. 2014; Bland-Hawthorn & Ger- 
hard 2016). 


3.6. Spherical Symmetry of SNRs 


The spherical symmetry of SNRs has been discussed 
by many authors in the past. Lopez et al. (2009) in 
their study from X-ray data presented evidence that 
SNRs can be be separated by SNe type (Ia or CC). 
Similarly, the infrared (IR) study done by Peters et al. 
(2013) shows that the SNRs resulting of a type Ia SN 
have a more circular and mirror symmetric morphol- 
ogy. However, as stated by (Dubner & Giacani 2015) 
and demonstrated by Ranasinghe & Leahy (2019), there 
is no discernible separation of SNR types according to 
their radio morphology. This may be caused hy the com- 
plex interactions between the SNRs and the ISM. 

We have examined whether there is a correlation be- 
tween the ovality and other SNR parameters. Here the 
ovality is estimated as defined by Bozzetto et al. (2017) 
(their equation 9). We found that none of the SNR pa- 
rameters (spectral indices, SNR ages, MC associations 
and galactic heights) or the explosion type was corre- 
lated to the spherical symmetry of SNRs. 


4. CONCLUSIONS 


We have compiled a table of 390 Galactic radio SNRs 
and their basic parameters. The following conclusions 
are formed from our analyses. 


1. The arithmetic mean of the Galactic SNR diame- 
30.5 pe with standard error 1.7 pe and stan- 


dard deviation 25.4 pe. The geometric mean and 
geometric standard deviation factor of the Galac- 
tic radio SNR diameters is 21.9 pe and 2.4, ге- 
spectively. While the arithmetic mean of Galactic 
SNR diameters is different from the mean SNR di- 
ameters from other galaxies, the geometric mean 
of Galactic SNR diameters is comparable to the 
SNR diameter means of M83. 


2. The E-D relation for the Galactic SNRs shows a 
large scatter in the X-D plane because the radio 
brightness is affected by a large number of physical 
factors (Pavlović et al. 2018). We did not find the 
scatter to reduce for particular classes of SNRs. A 
comparison of distances obtained by different X- 
D relations to distances obtained by more reliable 
methods show a discrepancy of 50% on average, 
questioning the reliability of the X-D method. 


3. We have estimated ages to 97 SNRs and on average 
D = (9.5240.92) (9-415901 pc for shell-type SNRs. 
However, this relation has a large scatter of about 
factor 2 (Figure 5). 


4. A birthrate of 1 SN in 58 yr fits the SNR age dis- 
tribution for ¢ < 2 kyr and is consistent with the 
currently accepted SN birth rate within 20. As 
age (t) increases, the mean birthrate drops, indi- 
cating the incompleteness of the sample of SNRs 
increases with age. 


5. The mean spectral index of the Galactic shell-type 
SNRs is —0.51 + 0.1. No correlations were found 
between the spectral index values and other SNR 
parameters (Le. MC associations, SN type, di- 
ameter and surface brightness). On average, the 
younger SNRs have a steeper spectral index. 


6. The majority of the SNRs are located near the 
Galactic plane and the Galactic height distribution 
of SNRs is best described by an exponential with 
scale height 48 + 4 pe. 


= 


'. The spherical symmetry of a SNRs in radio is not 
correlated to any other SNR parameter or the ex- 
plosion type. 


We thank the referee for the comments and sugges- 
tions that have improved this work, 

We acknowledge the support of the Natural Sciences and 
Engineering Research Council of Canada (NSERC). 
Nous remercions le Conseil de recherches en sciences na- 
turelles et en génie du Canada (CRSNG) de son soutien. 


п 


"ric SNRs 


A STATISTICAL ANALYSIS OF GALAC 


spev лора паэ aq st 922 pus spa sofas эца st 


irf HONDAS ur poquosp su saanreawd ратите Wm poymunsa ama. 


EMmIT-TLL 


is aya cere fouly = Gava чадра aesaat o - > 


cape sosy -q 


‘ausodino:) э р о роды “a Б sod UNS - vanos 
sr то ET Sweat Torri S 
" mo ero- т ахи o i0-rUD PE 
or dA zoo T E TxE OT-TTD dE 
oz ГТО En A o nl DIM ет вхи OO+OTID а 
o E DI 60 эхо OO+S0ID T 
o soo E PA РГ soror worse vo ro хат 30-662 OF 
o eros m аха 90-965 Gc 
le rore E эх о voit сотр  roFyo Le MXS 00260 sc 
o so 6 махра  rotesp uz 
o occu wA — Зоб кот Котев so- og sxs го 280 90 
D г vr S хор PS-LS cc 
te vor E о от Deer vt 5 pxo ое rz 
o SPFP di 026016: — rlcét Тор и $ ка KE-L ва 
o oe PIDE su S ахар rOIiD Gl 
o зготто- isz Б хат rO-0zD gc 
o im тож та Б віхар ro-ceD oc 
o ro- $ их op+ro GL 
D тазом оо se әд SOFTER LIII FOFSI 9г0= 90— о xst  rO-r9D sr 
ui [S я хо стъгео и 
or E 60- S шхяг cotTeD oF 
D iro- s — xoc e 
zz oo EZ SA дет Lore ror so oe Б ах PL 
o «тже 00 E SA КОРЕ то zo- а of er ет 
o so Er архаг 95-250 cr 
o [am t S хз T9910 TI 
o я sino Worsoe- 05 n вт S ехе soe OF 
o so- we $ ма Se- 6 
o roso ie Б ахар 0-915 9 
" S007 s90- vu S "XH TO-IfD L 
o oSI "I. dooT€rO0 A воза so oroms svo vo 5 roteto 9 
o 9A vs ES 5 ro-ro € 
o rozse E котом sro aroze ¿so sr s To-09 v 
o 2016 or PA voco sro oro H's L007 890- ist о sXs — rote £ 
19 ТЕТЯ E wo oro s гожо a 5 ахи 004500 c 
se rores хар zt PA o oro s oro- wor 5 crc ootoo т 
(oa) (аля) (oa) (od) d) (se) 
% sexu sug spei аку yuna sawap sous о хәрщ  Ausuop xng „dÂ 
Aymo  ampmg  aTwatay NS pay ow zip onus umop  pupadg m91 uns в us # 


siojoureged api pue SHNS WPO “T мачт, 


12 RANASINGHE & LEAHY. 


REFERENCES 


Anderson, L. D., Wang, Y., Bihr, S., et al. 2017, A&A, 605, 
A58, doi: 10.1051/0004-6361/201731019 
Badenes, C., Maoz, D., & Draine, B. T. 2010, MNRAS, 
407, 1301, doi: 10.1111/j.1365-2966.2010.17023.x 
Bland-Hawthorn, J., & Gerhard, О. 2016, АВАКА, 54, 
529, doi: 10.1146/annurev-astro-081915-023441 
Bozzetto, L. M., Filipović, M. D., Vukotić, B., et al. 2017, 
ApJS, 230, 2, doi: 10.3847/1538-4: 
Brogan, C. L., Gelfand, J. D., Gaensler, B. M., Kassim, 
N. E., & Lazio, T. J. W. 2006, ApJL, 639, L25, 
doi: 10.1086/501500 
Case, G. L., & Bhattacharya, D. 1 
doi: 10.1086/306089 
Chevalier, R. A. 1982, ApJ, 258, 790, doi: 10.1086/160126 
Cioffi, D. F., McKee, C. F., & Bertschinger, E. 1988, ApJ, 
34, 252, doi: 10.1086/166834 
Clark, D. H., & Caswell, J. L. 1976, MNRAS, 174, 267, 
doi: 10.1093/mnras/174.2.: 
Dopita, M. A., Blair, W. P., Long, K. S., et al. 2010, ApJ, 
710, 964, doi: 10.1088/0004-637X/710/2/964 
Dubner, G., & Giacani, E. 2015, A&A Rv, 23, 3, 
doi: 10.1007/500159-015-0083-5 
Ferrand, G., & Safi-Harb, S. 2012, Advances in Space 
Research, 49, 1313, doi: 10.1016/j.asr.2012.02.004 
Filipović, M. D., Payne, J. L., Reid, W., et al. 2005, 
MNRAS, 364, 217, doi: 10.1111 /j.1365-2966,2005.09554.x 
Green, D. A. 1991, PASP, 103, 209, doi: 10.1086/132810 
— 2005, Mem. Soe. Astron, Italiana, 76, 534. 
https://arsiv.org/abs/astro-ph/0505428 
—. 2019, Journal of Astrophysics and Astronomy, 40, 36, 
doi: 10.1007/s12036-019-9601-6 
Leahy, D. A. 2017, ApJ, 837, 36, 
doi: 10.3847/1538-4357/aa60el 
Leahy, D. A., Merrick, F., & Filipović, М. 2022, Universe, 
8, 653, doi: 10.3390/tniverse8120653 
Leahy, D. A., & Ranasinghe, S. 2018, ApJ, 866, 9, 
doi: 10.3847/1538-4357 /aade48 
Leahy, D. A., Ranasinghe, S., & Gelowitz, M. 2020, ApJS, 
248, 16, doi: 10.3847/1538-4365/ab8bd9 
Leahy, D. A., & Williams, J. E. 2017, AJ, 153, 
doi: 10.3847/1538-3881 /ааба 5 


/aa653c 


8, ApJ, 504, 761, 


239, 


Lee, J. H., & Lee, M. G. 2014, ApJ, 786, 130, 
doi: 10.1088/0004-637X /786/2/130. 

Long, К. S., Blair, W. P., Winkler, P. F., et al. 2010, ApJS, 
187, 495, doi: 10. 1088/0067-0049/187/2/495 

Lopez, L. A., Ramirez-Ruiz, E., Badenes, C., et al. 2009, 
ApJL, 706, L106, doi: 10.1088/0004-637X /706/1/L106. 

. Z., Dobardzié, A., Vukotić, B., & Urošević, D. 

2014, Serbian Astronomical Journal, 189, 25, 

doi: 10.2208/8A.J1489025P 

. Z., Urošević, D., Arbutina, B., et al. 2018, 
ApJ, 852, 84, doi: 10.3847/1538-4357 /ава1еб 

Peters, C. L., Lopez, L. A., Ramirez-Ruiz, E., Stassun, 
К. G., & Figueroa-Feli 2013, ApJL, 771, 
doi: 10.1088/2041-8205/771/2/L38 

Ranasinghe, S., & Leahy, D. 2019, arXiv e-prints, 
arXiv:1909.11803. https://arxiv.org/abs/1909.11803 

—. 2022, ApJ, 940, 63, doi: 10.3847/1538-4357 /ac940a. 

Roberts, М. Е. 2004, The Pulsar Wind Nebula Catalog 
(March 2005 version), 
https:/ /www.physies.megill.ca/~pulsar/pwneat. html 

Rozwadowska, K., Vissani, F., & Cappellaro, E. 2021, 
NewA, 83, 101498, doi: 10.1016/j. 20.101498 

Sedov, L. I. 1959, Similarity and Dimensional Methods in 
Mechanies 

Бирав, L., Castelletti, G., Peters, W. M., & Kassim, N. E. 
2018, A&A, 616, A98, doi: 10.1051 /0004-6361 /201832995 

Tammann, G. A., Loeffler, W., & Schroeder, A. 1994, 
ApJS, 92, 487, doi: 10.1086/192002 

Urquhart, J. S., Figura, C. C., Moore, T. J. T., et al. 2014, 
MNRAS, 437, 1791, doi: 10.1093 /mnras/stt2006 

Wang, Q. D., Lu, F. J., & Gotthelf, E. V. 2006, MNRAS, 
367, 937, doi: 10.1111/j 2966.2006.09998.x 

Weiler, K. W., & Panagia, N. 1978, A&A, 70, 419 

Winkler, P. F., Coffin, S. C., Blair, W. P., Long, K. S., & 
Kuntz, K. D. 2021, ApJ, 908, 80, 
doi: 10.3847/1538-4357/abd77d 

X-Z, & Han, J.-L. 2005, ChJA&A, 5, 

165, doi: 10.1088/1009-9271/5/2/007 


Pavlović, 


Pavlov. 


8, 


wast. 20 


